Effective three-particle interactions in low-energy models for multiband systems 
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We discuss different approximations for effective low-energy interactions in multiband models for 
weakly correlated electrons. In the study of Fermi surface instabilities of the conduction band(s), 
the standard approximation consists only keeping those terms in the bare interactions that couple 
only to the conduction band(s), while corrections due to virtual excitations into bands away from 
the Fermi surface are typically neglected. Here, using a functional renormalization group approach, 
we present an improved truncation for the treatment of the effective interactions in the conduction 
band that keeps track of the generated three-particle interactions (six-point term) and hence allows 
one to include important aspects of these virtual interband excitations. Within a simplified two- 
patch treatment of the conduction band, we demonstrate that these corrections can have a rather 
strong effect in parts of the phase diagram by changing the critical scales for various orderings and 
the phase boundaries. 

PACS numbers: 74.20.-z, 74.62.-c 
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I. INTRODUCTION 

In various areas of recent condensed matter physics it 
has become clear that interaction effects in multiband 
systems can lead to interesting phenomena that usually 
do not occur in single-band models. An important ex- 
ample thereof is the gap structure of the iron pnictide 
superconductors, where calculations show that the gap 
structure of the superconducting state depends signifi- 
cantly on the intra- and interorbital interaction param- 
eters, and on the orbital composition of the bands that 
are usually mixtures of several Fe d-orbitals (see, e.g., 
Refs. [l|-l5|)- If this complexity was ignored, one would 
obtain rather isotropic gaps around the Fermi surfaces. 
The gap anisotropy arises due to orbital makeup^, i.e. 
due to the additional wavevector dependence of the or- 
bital components of the Bloch functions for the bands 
that give much more structure even to local interactions 
in these multi-orbital systems. 

The additional structure in the multiband interactions 
also creates a stronger sensitivity to features and param- 
eters of the respective effective low-energy model, com- 
pared to one-band models. This calls for a more compre- 
hensive theoretical study of multi-orbital model systems. 
One important aspect is whether the orbital composi- 
tion of the bands is altered by interaction or correlation 
effects beyond those already captured in density func- 
tional theory (DFT). Another aspect is that all these 
studies are performed within effective low-energy models 
where many bands outside a respective energy window 
have been absorbed into model parameters. The ques- 
tion that shall be in the focus of this work is if there are 
sizable additional corrections due to virtual excitations 
in the bands outside the low-energy window in which the 
effective model is formulated. In particular we will try 
to estimate if these virtual excitation effects can be more 
important than orbital makeup effects. 

A concrete physical question that motivates such con- 



siderations occurs in the high-T c cuprates. Here a com- 
parison of DFT electronic structure data and experimen- 
tally measured critical temperatures T c for superconduc- 
tivity suggests^ that lowering a so-called axial orbital in 
energy from above toward the Fermi-surface-forming Cu 
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-orbital causes an increase of T c , although the 



same change makes the Fermi surface more round and 
hence less favorable for antiferromagnetic spin fluctua- 
tions. In a recent studyi we have tried to assess if orbital 
makeup effects can cause such material trends, but found 
only very moderate improvement of T c . Hence the ques- 
tions is if the next step in a reduction of approximations 
results additional increases of T c . 

We already note at this stage that virtual-excitation 
corrections have been identified as an important mech- 
anism to screen down local (Hubbard) interaction pa- 
rameters. A currently popular scheme to capture this 
effect is the constrained-RPA approach^ that sums up 
the RPA series for the Coulomb interaction with at least 
one intermediate line in the bubble being a high-energy 
excitation. In another way, the work described here can 
be understood as an attempt to partly include such con- 
tributions in a renormalization group approach for the 
effective low-energy model. 

In order to detect superconducting instabilities and 
hence to estimate a superconducting energy scale the- 
oretically, one should look for a divergence of the effec- 
tive interactions at low scales. In one-loop RG studies 
this is accomplished by integrating out the single-particle 
modes of the effective model with a decreasing cutoff or 
RG scale (for a recent review, see Ref . [lot) . In this proce- 
dure, one-loop corrections with both internal lines in the 
effective low-energy window are summed up. The initial 
two-particle interaction for this flow should include the 
influence of the excitations outside the low-energy win- 
dow. We can try to estimate these corrections in second 
order perturbation theory, which gives rise to diagrams 
with two propagator lines connecting two interaction ver- 
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tices. For these diagrams we will have contributions with 
the two internal lines both in the high-energy range, and 
with one in the high-energy and the other in the low- 
energy window. The latter processes can be expected to 
be more important, as they usually come with the smaller 
energy denominator. In the following we will argue that 
systematically including these latter contributions into 
the solution of the effective low-energy problem requires 
either an approximate correction of the low-energy effec- 
tive four-point vertex (two-particle interaction) or to im- 
prove the truncation so as to keep the effective six-point 
vertex (three-particle interaction) in the effective low- 
energy model. We will compare these two schemes quan- 
titatively in a simplified model and demonstrate that the 
virtual excitation effects can indeed change critical scales 
for superconducting instabilities considerably. 



Upper band 




FIG. 1. (Color online) Band structure of the two-band model 
for AE — 19.5, Wd = 11.5 and w s = 4.5 at van-Hove filling. 
The lower band on the right is the conduction band. The 
colorbars on the sides indicate the band energy. This choice 
of parameters underlies the RG flow in Fig. [6] 



II. MODEL 

The basic model for this work is a two-orbital Hamil- 
tonian where one of the two resulting band crosses the 
Fermi level, while the other one is separated by an energy 
gap. We will discuss several approximate ways to include 
the band away from the Fermi surface into effective one- 
band models for the band near the Fermi surface. The 
structure of the model we use is borrowed from an (ef- 
fective) two-orbital model derived for high-T c cuprates, 
more precisely for the spa and dpa orbitals from the 
local density approximation (LDA) band structure of 
YBa 2 Cu30?ii, extended by inter- and intra-orbital in- 
teractions of strengths U and U', 
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Here, /k,a,er and n iya ^ = fl a ^fi., a ,a denote the annihila- 
tion operator of an electron with momentum k, orbital a 
and spin orientation a and the density of such a fermion 
at site i, respectively. The non-interacting part of the 
Hamiltonian is given by 

A k = AE + (1 - u+) Ws , 
Bk = (l - u^) Wd , 
Ck = -u k y/w d w s , 

AE being the band separation and w s and Wd the widths 
of the s- and d-orbitals, respectively. Of course other in- 
teractions terms like a Hund- or non-local terms can be 
added, but this does not play any role for the considera- 
tions that follow. 



Our motivation for studying this model is that serves 
as a simple test case for the development of an func- 
tional renormalization group (fRG) approach to multi- 
band systems. We do not aim at making predictions for 
a specific material that hold on a quantitative level. Re- 
garding YBa2Cu307, our results should only be taken 
with a grain of salt. The reason for this is twofold: 
First of all, a multi-orbital tight-binding model of the 
cuprates should include orbitals on the oxygen atoms. 
The two-orbital model given above, however, does not 
allow for the description of Varma currents^ and other 
types of intra-unit cell order—. Moreover, the fRG ap- 
proach for fermions used in this work is only viable as 
long as the renormalized interaction stays weak, whereas, 
in the cuprates, realistic values for the bare interaction 
are already large compared to the bandwidth. Also other 
parameters will not always be chosen according to ab ini- 
tio calculations in the following. 

We now switch to an imaginary-time functional inte- 
gral formalism with Grassmann fields ipk,a.a and ipk,a,a- 
corresponding to the operators a and fu,a,a- These 
fields depend on the 1 + 2 momentum k — (fco,k) with 
Matsubara frequency ko. 

Diagonalization of the quadratic part of the action 
corresponding to Eq. (|TJ) yields bands with energies 
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which are shown in Fig. Q] We will work in a parameter 
range where the lower of the two bands cuts the Fermi 
level, while the upper band is entirely above the Fermi 
level. The effective low-energy model will the refer to a 
model with degrees of freedom in the functional integral 
that reside in the lower band near the Fermi surface only. 
This simple two-band set-up serves as minimal model for 
more complex situations in with many, possibly entan- 
gled bands in both high-energy and low-energy sector. 

The equation for the unitary transformation from or- 



bitals i[) to bands \ labeled by a = ± reads as 

1pk,a,tr = ad k Xk,a,cr + C k Xk,~a,a , (2) 



Ck = N k Ck, 

where iVk normalizes the transformation to a unitary one. 
The inverse of this transformation gives the orbital am- 
plitudes adk and c k for the band fields Xk,a,a- If we ex- 
press the local interactions that are quartic in the fermion 
fields in terms of the band fields using Eq. @, we obtain 
four factors of dk, a ,cr or Ck, a .a multiplying the interac- 
tion parameters U and U'. This collection of prefactors 
was dubbed orbital makeup! and obviously leads to some 
wavevector-dependence even of local interactions when 
they are expressed in band language. In Appendix lAl we 
parameterize the interaction in the band picture in an 
efficient way, and later we will also discuss the impact of 
this orbital makeup on the critical scales for supercon- 
ducting instabilities. 



III. EFFECTIVE LOW-ENERGY ACTION 

Now we proceed towards an effective one-band action 
for the band near the Fermi level. This means we want to 
integrate out the upper band £7k,+ away from the Fermi 
level. In order to be concise, we use the matrix notation 

4>Mr] = ^2 <f>k,a,(TMk,a,*;k> ,a> ,<r'»?fc' ,a> ,<r> 

fc,a,cr k' ,a' ,ct' 

for the free part and other field-bilinears of the action. 

Consider now the generating functional W of the con- 
nected Green's functions 



W[q. rj\ 



In 



-S[x,x] e vx+xv 



with source fields 77. We are interested in the low-energy 
properties of the system. Before taking derivatives with 
respect to the source fields, their upper-band components 
can thus be set to zero. In the band language, the fields 
and the quadratic part S^ 2 ' of the action can be decom- 
posed into lower and upper band parts 



S {2) [X,X] = XDX = X+D+X+ + X D 



with 



D = D A 
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Therefore, the covariance splitting formula 14 - applies, 
which gives rise to the following form of the one-band 
effective action S e g for the lower band 

S e $[X-,X-} = X D x +V[x"_,x_] , 

\ J = lj) x+e -X+D + X +e -S W [x++X-,X++X-] ^ 

(3) 








FIG. 2. Four- and six-point vertex of S e ff for the effective 
one-band model as considered at approximation level 2. Small 
filled vertices correspond to the bare interaction S^ 4 ' , whereas 
the vertices of S e a on the left hand sides are denoted by empty 
circles. Solid lines represent the lower and dashed lines the 
upper band. There are no propagators attached to the exter- 
nal legs. Self-energy effects will be neglected. 



which leads to 
W[fj-,T]-] = - In 



T>x- e 



-S e{ f[x-,X-] e V-X-+X-V- 



The effective interaction V contains a functional integral 
over the upper band part of the fields with measure T>x+ 
and corresponds to the generating functional of ampu- 
tated connected Greens functions, as used in the Polchin- 
ski renormalization group scheme^ (for a comprehensive 
reviews of the various generating functionals in our con- 
text, see e.g. Refs. Il6l - [l8l ). This means that the parame- 
ters of the effective low-energy action are given by these 
amputated connected Green's functions. In the special 
case of Eq. however, only the upper band has been 
integrated out. Thus, in the diagrammatic expansion of 
the expansion coefficients in the fields X- and x_, the 
propagators on internal lines are restricted to the upper 
band, whereas external legs live on the lower one. Before 
we evaluate S e ff, we recall that the amputated connected 
Greens functions can be recovered from one-particle irre- 
ducible (1PI) diagrams by drawing all tree diagrams with 
1PI vertices. In our case, the internal lines of these tree 
diagrams are high-energy propagators. 

For assessing the low energy properties of the two-band 
model, the general strategy is as follows. In a first step, 
we calculate S c f{. Since the Fermi surface does not inter- 
sect with the upper band, the diagrams that appear in 
Eq. © need not be regularized and therefore the lower- 
band effective action can be evaluated perturbatively for 
sufficiently small values of the bare interaction. In a sec- 
ond step, the effective action of the lower band is treated 
by a method of choice, which is in our case the fRG 10 . 
Then S c s imposes an initial condition on the RG flow. 

For our main focus, determining the energy scale for 
Cooper instabilities, the effective interaction of the low- 
energy model is the object of prime interest. The sim- 
plest truncation of the effective action would then be 
to drop all terms of order higher than four. Then the 
four-point term in the effective action, S^S , can be com- 
puted in perturbation theory in the interactions involving 
high-energy modes. In lowest order, i.e. zeroth order in 
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the interactions with high-energy modes or first order in 
the bare interactions irrespective of energy scales, S^a is 
just the bare interaction of the fields in the lower band, 
decorated with orbital makeup. This is the standard 
that has been used e.g in all studies of unconventional 
pairing in the iron pnictides, or is used implicitly when 
the one-band Hubbard model is used as a model for the 
high-Tc cuprates. We will later denote this truncation 
as approximation level 1. Solving the low-energy theory 
diagrammatically captures (possibly singular) diagrams 
with both internal lines in the low-energy window. 

In the next order for , i.e. second order in the bare 
high-energy interactions, we get various diagrams. On 
one hand, there are self-energy Hartree- and Fock-like 
contributions on the external legs. Further, there are one- 
loop corrections with both lines in the high-energy sector. 
As mentioned, these are non-singular one-loop terms, as 
all internal lines are away from the Fermi level. Let us 
call the scheme that keeps these terms approximation 1 '. 
If we now again solve the low-energy theory diagrammat- 
ically, we will capture one-loop corrections for the effec- 
tive interaction that have both propagator lines either in 
the high-energy range and that are already included in 
, and corrections with both lines in the low-energy 

sector, coming from the perturbation expansion in S^g ■ 
What is not included are 'mixed' diagrams with one in- 
ternal line in the high-energy range and one line in the 
low-energy window. Looking at the energy denominators 
of these lines, these excluded mixed contributions should 
be potentially more important than those with two inter- 
nal high-energy propagators captured at approximation 
level 1'. 

We could now go on and include further order in the 
bare interactions as corrections to Sj$ . It is however 
clear that in these higher-order diagrams all internal lines 
will be high-energy propagators. Hence, these corrections 
do not include the missing mixed diagrams, and thus 
we do not pursue these corrections any further here. In 
principle, they can be summed up using RG schemes, as 
described in Ref. [H. 

The next useful extension would be to go along an- 
other path and to improve the truncation of the effective 
action. Hence we now keep the six-point term in the ef- 
fective action. In the tree diagram expansion of S s the 
sixth-order term is generated in second order in the 1PI 
four-point vertices of the high-energy theory. Here we 
will replace them by the lowest order, i.e. by the bare in- 
teractions, as shown in Fig. [3] This means that possible 
renormalizations of the two-particle interactions by ad- 
ditional high-energy processes are deliberately excluded. 
As argued above these corrections with additional high- 
energy propagators should however be smaller due to the 
energy separation of the bands. This we will call approx- 
imation level 2. Furthermore, we should potentially drop 
self-energy corrections in order to avoid double counting 
of contributions, i. e. to the self-energy that are already 
included the DFT-based derivation in our two-orbital 



Hamiltonian. In this approximation, the quadratic part 
and the bare four-point couplings remain unrenormal- 
ized, whereas a six-point term depicted in Fig. (J2|) is 
generated. For the effective action, we get 

S cS [X-,X-} = (S (2) + S (4) ) {X-,X-} +S^{x-,X-} 

xx-(&)x-(&)x-(&), 

d = (ki,CTi) being a short-hand notation for the quan- 
tum numbers of the fields. For the precise form of the 
six-point coupling function we refer to Eq. (|B1[) in 
the appendix. Now, in a perturbative expansion for the 
low-energy theory, we will receive contributions where 
two legs of the six-point coupling will be folded together 
by a low-energy propagator line. As the six-point term 
came about by joining two legs of two four-point inter- 
actions by a high-energy propagator line, this will effec- 
tively bring in those missing diagrams with two internal 
propagators, one of which is a high-energy mode, and the 
other a low-energy mode. 

We note in passing that the constrained RPA used for 
computing effective Hubbard interaction parameters^ 
can be understood as an infinite order resummation of 
the mixed diagrams included at approximation level 2. 
Resummation at level 1', i.e. without any internal lines in 
the low-energy sector, would presumably result in much 
less reduction of the onsite repulsion. On the other hand 
we note that just keeping the six-point term still does not 
capture the full cRPA series, as pure powers of mixed 
loops (i.e. bubbles with one high-energy and one low- 
energy line) included in the cRPA are not contained in 
the RPA series generated at our approximation level 2. 
This can be seen from constructing bubble sums with 
the elements of Fig. [2] by contracting low-energy lines. 
There is a mixed diagram in second order in the bare 
interactions, but in third order or fourth order we have 
to add a pure low-energy loop between two mixed loops. 
So, some orders in the mixed diagrams are still missing, 
but including the six-point term goes in the right direc- 
tion. On the other hand, in contrast with the cRPA, 
our approximation does not neglect vertex corrections or 
particle-particle diagrams. The goal of the present study 
is to show the effects of these corrections to the simpler 
truncation. The question of how the cRPA series is un- 
derstood in terms of the effective action will be discussed 
in another forthcoming publication. 



IV. FRG TREATMENT 

We are now in a position to proceed with the second 
step of solving the low-energy model, which we will do 
by a fRG flow. This will clarify the differences between 
the various levels of approximations. 
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FIG. 3. Flow equations for the self-energy, the four-point and 
six-point vertices, all one-particle-irreducible. A dot repre- 
sents a derivative with respect to the cutoff. Lines with a slash 
correspond to the single-scale propagator S = G — GSG. For 
more details, see e.g. Ref. ITol . 



A. 1PI functional RG scheme with smooth 
frequency cutoff 

We now introduce an infrared (IR) cutoff on the lower 
band by replacing £>_ by D^R^ 1 , where R\ denotes a 
regulator function. In this paper, we choose 

This particular choice of the regulator, introduced by 
Husemann and Salmhofer 20 , does not completely sup- 
press contributions from the Fermi surface at nonzero A 
and therefore allows one to take a possible ferromagnetic 
instability into account. Moreover, a pure frequency cut- 
off with R\ = at fco = circumvents Fermi-surface- 
renormalization issues^, since the full propagator reads 
as 

G = R X [Q + HRaF 1 , 

X being the self-energy. So self-energy effects may be 
neglected without ignoring the most relevant terms. 

Starting from an exact flow equation for the generating 
functional T_ of the one-particle irreducible (1PI) vertex 
functions related to W[fj_ , ry_] by a Legendre transforma- 
tion in the lower-band fields, one obtains an infinite hier- 
archy of differential equations for the 1PI vertice o 10 ' 14 . 
The subscript of T_ reminds us that only degrees of 
freedom in the lower band are integrated out, while the 
higher bands have been absorbed in the initial conditions. 
The first three of the RG equations for the vertices gen- 
erated by r_ are expressed in a diagrammatic language 



in Fig. |3l At this point we note that there is no sim- 
ple relation between T_ and the Legendre transform T of 
W [77, 77] . In the former case, information about correla- 
tions in the upper band is lost and the one-particle irre- 
ducibility only holds regarding propagators on the lower 
band. In the latter, namely the full multiband case, the 
band index however appears as an additional quantum 
number that is summed over on the internal lines of f PI 
diagrams. 

In order to make progress, the hierarchy of flow equa- 
tions needs to be truncated at some point. In the con- 
ventional truncation scheme one neglects the six-point 
vertex completelyi 4 -. In the so-called Katanin schemed 
six-point contributions that are generated during the RG 
flow are fed back into the flow equation for the four-point 
vertex. Both established truncation schemes however are 
not suited for a non- vanishing initial six-point vertex. So 
its impact on the flow poses a conceptually new problem. 



B. Improved truncation with one- loop six-point 
feedback 

Let us now return to the initial interactions of the low- 
energy problem in the conduction band that are given 
by the effective interactions after the high-energy modes 
in the upper band have been integrated out. As argued 
above, dropping all effective interactions higher than the 
four-point (two-particle) term ignores possibly important 
contributions. Keeping the six-point term in arbitrary or- 
der in the bare interactions or even keeping higher terms 
are hard tasks. A simpler, improved truncation scheme 
going beyond approximation level 1 consists in neglect- 
ing eight-point and higher interactions and assuming that 
the six-point vertex does not flow in the weak coupling 
regime. Then the six-point vertex would just be given by 
the product of two 1PI four-point vertices, connected by 
a high-energy propagator as depicted in the diagram in 
the lower part of Fig. [2j If the six-point vertex is then fed 
back into the four-point flow equation for the low-energy 
theory, the missing diagrams appear on the right-hand 
side with one high-energy and one single-scale low-energy 
line as shown in Fig. |4j In this truncation, two-loop con- 
tributions are neglected. Note that in Ref. also two- 
loop terms have been considered. However, this has only 
been done for the case of an initially vanishing six-point 
vertex and purely local bare interactions. 

Before solving the flow equation in this truncation let 
us look briefly at the feedback term consisting of the first 
two diagrams in Fig. 2J Their precise form shall be given 
in Appendix [Bj The fact that the diagrams with self- 
energy insertions in Fig. U get disconnected when the line 
on the upper band is cut should not lead to confusion. 
Since we are dealing with an effective lower-band action, 
one-particle irreducibility only holds regarding lines on 
the lower band. If self-energy effects are neglected, the 
extra term in the flow equation for the four-point ver- 
tex is just identical to the scale derivative of the sum of 




FIG. 4. Effective ffow equation in the low-energy problem 
for the one-particle irreducible four-point vertex including the 
feedback term re-expressed in terms of the bare couplings. 
The dashed line stands for propagators in the upper band. 



second-order diagrams in the bare interaction that have 
one internal line one the upper band and another on the 
lower one. In the infrared, they stay regular but will be 
larger compared to second order diagrams with all inter- 
nal lines on the upper band provided that the band sep- 
aration is sufficiently large. We therefore neglect those 
upper band diagrams and restrict the terms in the ef- 
fective interactions to the tree level in the high-energy 
modes. 

In the following, we shall distinguish three levels of 
approximation: 

1. Ignoring the six-point vertex completely. In this 
approximation the only multiband effects are the 
signatures of orbital makeup. We consider this con- 
ventional truncation for comparison in order to dis- 
tinguish orbital makeup from six-point effects. 

2. Including the feedback term in one-loop fRG and 
using the flow equation Fig. H) 

3. Adding the mixed-band diagrams in the limit A — > 
to the initial condition for the flow of the low- 
energy model. This approximation should yield 
reliable results if the mixed diagrams are already 
close to their infrared value at scales at which the 
lower-band diagrams only have induced a small 
renormalization of the initial couplings. 

Approximation levels 1 and 2 were already introduced 
in Sec. IIII1 and level 3 is a simplification of level 2 that 
is easier to handle numerically. In the following we will 
compare the fRG flows in the low-energy models resulting 
from these approximations. 

If one recalls that the LDA-derived dispersion of the 
underlying two-orbital model already contains interac- 
tion effects on a certain level, the band-flip self-energy 
insertion diagram (second term on the right hand side in 
the diagrammatic equation in Fig. [4]) should potentially 
be neglected in order to avoid double counting. Its im- 
pact is however not important, as will be commented on 
at the end of Subsec. IV Bl 
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FIG. 5. The four running couplings g\ to <?4 of the two-patch 
approximation for the conduction band. The blue lines rep- 
resent the Fermi surface. 

C. Two-patch approximation 

In order to make the resulting numerical calculations 
more feasible, and as we are mainly interested in getting a 
first picture of the effects due to the higher truncation, we 
now employ the so called two-patch approximation that 
has been used in the context of the one-band Hubbard 
mode l 23 i 24 or iron pnictides^. 

Just like the one-band Hubbard model the dispersion 
of the lower band of our model Hamiltonian in Eq. ([T]) 
has saddle points at k = A = (0,7r) and B = (it, 0). If 
the system is now considered at van Hove filling, A and 
B lie on the Fermi surface. At zero temperature, the 
low energy properties then are dominated by contribu- 
tions of a small vicinity around these saddle points. We 
therefore restrict the internal momenta in the lower band 
diagrams to two small patches around A and B. In the 
following, we neglect self-energy effects and the frequency 
dependence of the coupling functions. In the one-loop di- 
agrams the momentum integral then only enters in the 
bare susceptibilities 

*pp(0 = E / dpG-(p)G-(l- P ) 

Po "'patch 

<M0 = E / dpG-(p)G-(l+p). 

p "'patch 

For zero temperature, the Matsubara sum is evaluated 
analytically for Iq = and the two-patch approximation 
restricts the transfer momenta 1 to or 7r = (7r, it). In 
our truncation, the fRG analysis can now be restricted 
to four running couplings depicted in Fig. [5j namely 

gi = V-(A,B,B) = V-(B,A,A) 
g 2 = V- (A, B,A) = V- (B, A, B) 
g 3 = V-(A,A,B)=V-(B,B,A) 
9i = V-(A,A,A) = V-(B,B,B), 

which drastically reduces the computational cost of the 
RG flow. At approximation level 1, the initial conditions 
for these couplings are obtained from transforming the 
intra- and interorbital interactions from the bare Hamil- 
tonian into the band representation. The corresponding 
equations are given in Appendix IA1 Eq. (|A2[) . 

As for the mixed-band diagrams, however, also re- 
gions away from the Fermi surface contribute signifi- 
cantly. Therefore, the loop integrals in the mixed-band 
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diagrams are to be taken over the whole Brillouin zone. 
In the case of the one-band Hubbard model, the cut- 
off can be chosen such that the resulting loop integrals 
can be evaluated analytically and do not depend on the 
patch size2£. However, that cutoff scheme is only vi- 
able in a small neighborhood around the saddle points 
whereas the cutoff needs to be defined on the entire Bril- 
louin zone in order to consider the six-point feedback. 
To the authors' knowledge, only momentum shell cutoff 
schemes have been used for the two-patch model while a 
frequency cutoff is used in this paper. 

The flow equations in the two-patch approximation 
read as 

<h = SV- (A, B, B) + di (<?? + si) + 2d 2 (g 2 - 9l ) 9i 
-d 3 (gl+gl) (4) 

g 2 = SV- (A, B, A) + d l92 ( 9l - g 2 ) + 2d 2 g 29i 

- 2d 3 g l92 (5) 

.93 - SV- (A, A, B) - 2d gm + 2di3s (2 9l - g 2 ) (6) 

g 4 =SV-(A,A,A)-d (gl+ 9 l) 

+ d 2 (gl + 2g 1 g 2 -2g 2 1 +gl) , (7) 

where the dot denotes a derivative with respect to A. 
The six-point feedback leads to correction terms SV- that 
do not occur in two-patch studies of one-band systems. 
They are given in appendix [B] The integration over the 
patches in the loops 

d a = *p P (0) , d x = $ p h(7r) , 

d 2 = % h (0), d 3 = -4 pp (7r) 

are performed numerically using an adaptive routine 27 . 
The different levels of approximation introduced in Sub- 
sec. HVBI now imply the following: 

1. Neglecting SV- and initializing the giS by the re- 
spective values of the coupling function V- at A — 
oo. 

2. Keeping SV- and initializing the giS by the respec- 
tive values of V- . 

3. Neglecting SV- and initializing the giS by the 
respective values of V_ + AV-, where AV- = 
— J °° dX SV- denotes the sum of all second-order 
mixed-band diagrams. 

In all the examples considered in this paper, we observe 
an abrupt flow to strong coupling at some critical scale 
A cr it • In the two-patch approximation, the flow equations 
of the couplings to external source fields for s- and d-wave 
superconductivity (a s sc an d aysc, respectively), anti- 
fcrromagnetism («af) an d ferromagnetism (ofm) take 
the simple forms 

a s sc = ~2d (.93 + gi) a s sc 
«dsc = -2d (.94 - gz) adsc 

"AF = +2rfl (.91 + 33) OAF 

c*fm = +2d 2 (.92 + g^) qfm • 



These couplings also appear in the expressions for the 
corresponding susceptibilities and determine which or- 
dering tendency grows fastest at the critical scale, i.e. is 
the leading instability at A CI it . 



V. NUMERICAL RESULTS 

We now solve the two-patch model at the three approx- 
imation levels described before for the parameters given 
in Tab. U at van Hove filling. We first discuss the approx- 
imation level 1, i.e. without six-point feedback and then 
continue to describe the changes in the higher approxi- 
mation levels. 



A. Flows without six-point term: Impact of orbital 
makeup and competition of FM and dSC instabilities 

In this subsection we discuss the impact of orbital 
makeup when the six-point feedback is neglected. 

Let us first consider the case of vanishing inter-orbital 
interaction, U' = 0. In that case, according to Eq. (|A2[) 
all four couplings take on the same value U e s at the be- 
ginning of the flow. Since d\ + c 2 = 1, these couplings 
are smaller than the intra-orbital coupling U . Moreover, 
upon an expansion of the lower band dispersion around 
the saddle points, the lower-band dispersion is equivalent 
to the dispersion of the one-band t-t' Hubbard model up 
to second order with effective parameters t and t' for 
next nearest and next-to-nearest neighbor hopping. So 
for U' = and approximation level 1, where the up- 
per band does not enter, we are back to the one-band 
t-t' model in a two-patch approximation, albeit with a 
reduced U and employing a smooth frequency cutoff un- 
like the the momentum-shell cutoff in previous studies of 
this model. 

First we characterize the nature of the flow in the ab- 
sence of a six-point term. Characteristic curves for the 
flow of the couplings gi for this approximation level 1 
can be seen in Figs. |5] and [7] For zero next-to-nearest 
neighbor hopping, the Fermi surface is perfectly nested, 
giving rise to antiferromagnetism as the leading instabil- 
ity. If t' is small but nonzero, the most important loops 
in the flow equations Eqs. (HUT]) are do and d\. Therefore 
d 2 and d 3 can be neglected at first. In this approxima- 
tion (?4 is decreased by do while the d\ term in Eq. ([6]) 
prevents g 3 from being renormalized to zero. This al- 
lows for a sign change of 54 after which the growth of 
do drives the system to a d-wave superconducting dSC 
instability corresponding to a strong coupling fixed point 
with <?4 — > — 00 and 91,32,33 — * +00. This instabil- 
ity is in full agreement with the findings of previous 
two-patch studies^,. If we now increase t' , the particle- 
hole diagrams with small wavevector transfer become 
more important, i.e. the d2 terms cannot be neglected 
any more. In particular, the last term in Eq. ([7]) ham- 
pers the sign change of 34, which now occurs at a lower 
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scale, leading to a lower critical scale. Moreover, g 2 is 
now renormalized to zero instead of diverging to +00. 
Thus we have an altered strong-coupling fixed point, but 
still with 53 — » +00, 54 — > —00, corresponding to a 
(i-wave pairing instability. This distinction was absent 
in traditional two-patch studie o 23 ' 24 in which the small- 
wavevector-transfer particle-hole channel was disadvan- 
taged by the choice of the cutoff function. Katanin and 
Kampf have included such contributions in a momentum- 
shell approach 26 . They find a similar strong coupling 
fixed point with non-diverging <? 2 which however corre- 
sponds to an antiferromagnetic instability. In A-patch 
flows we will a smooth crossover from the fixed point 
with (72 — > 00 to the one with g 2 — > 0. 

If the dispersion is varied further by increasing t' , d 2 
grows even more strongly. Then it eventually prevents 
a sign change of g± and therefore excludes d-wave super- 
conductivity. Then the flow corresponds to a FM insta- 
bility. Since the d 2 terms in Eqs. (jU[5]) depend linearly 
on 54, its sign change to negative values occurring in the 
(i-wave regime prohibited ferromagnetism, which reflects 
the mutual exclusion of the FM and dSC instabilities. 
Along the separatrix between these two regimes, all four 
running couplings flow to zero. This suggests that the 
critical scale drops to zero from both sides, which implies 
the existence of a quantum critical line between the two 
phases. In an A-patch study, however, the situation is 
more involved and both instabilities may occur simulta- 
neously. Unfortunately, the region around the separatrix 
is unaccessible in our calculations due to an excessive 
number of function calls required for numerical integra- 
tion of the loops. In Fig. [5J we plot the phase diagram of 
this model obtained with approximation level 1. In the 
e?-wave regime a reduction of the interaction strength to 
Ucs leads to a lower critical scale. 

Let us now turn to the case of non-vanishing inter- 
orbital coupling, where we have g\ — #4 > g% — 53 at the 
initial scale. From Table [TT| we find that this detuning 
leads to a lower critical scale in the dSC regime, as the 
inter-orbital interaction suppresses 53 — g± in the initial 
condition while g% + g± remains unchanged. 

Before we analyze the ci-wave regime in further detail, 
we briefly have a look at the band structure for parame- 
ters given for YBa2Cu3C>7 in Ref. Ill For such a system, 
however, van Hove filling is not close to the experimental 
situation, since it corresponds to a filling factor of about 
0.19. The result should therefore not be taken as a realis- 
tic prediction for this material. In Fig. [7] we observe that 
feedback the flow approaches a FM fixed point at level 1. 
This also holds for the other approximation levels. 



B. Inclusion of the six-point term 

We now turn our attention to the six-point feedback for 
the band structure underlying the flows in Figs.l9land[T0l 
employing the improved approximation levels 2 and 3 
introduced in Sec. IIVBI Note that the model parameters 
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TABLE I. Different parameter sets for the dispersion 
analyzed in this paper. The lower-band filling factor 
J BZ rfk 0[— E- (k)]/(47r 2 ) is denoted by ni ow and the maximal 
hybridization matrix element by c max = niax^ cu . 
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TABLE II. Critical scale A cr it and leading instability for dif- 
ferent levels of approximation. 



will be varied in order to clarify the differences between 
these approximations and may not always correspond to 
an experimentally realistic situation. 

We observe a flow to a dSC strong-coupling fixed point 
over a wide parameter range for all approximation levels 
considered. This pairing instability occurs irrespective 
of the presence of inter-orbital interactions. The crit- 
ical scale, however, is enhanced by the six-point feed- 
back. Some numbers can be inferred from Table HTl and 
also Fig. [TTJ the enhancement can easily be a factor of 
two, at least in this parameter range of small instability 
scales. We observe that A cr it only differs weakly between 
the approximation levels 2 and 3. Figs. [9] and [10] illus- 
trate that first the mixed-band diagrams flow to a value 
close to their infrared limit before the lower-band dia- 
grams start to grow significantly. (All other ratios of the 
couplings that are not shown in Figs. [§] and [TU] behave 
indeed likewise.) This separation of scales is enhanced 
or might be even induced by a two-patch approximation. 
It ensures that the two-loop correction term discussed 
in Appendix [C] remains negligible and also holds in the 
ferromagnetic case and for Fig. [TT] 

The question now is whether the six-point feedback on 
the critical scale can be related to characteristic proper- 
ties of the band structure such as band curvatures. For 
this purpose, we consider AV_ terms in the initial condi- 
tion of approximation level 3 for small hybridization Cfc 
as in the case of Fig. If for simplicity U' is then sent 
to zero, the bare coupling functions V- and V3 Eqs. (|A2[) 
and (|A3[) read in leading order in the hybridization 

V-(k u k2,k 3 )nUY[d ki 

i 

V 3 (ki,k 2 ,k 3 ) w -Uc kl d k2 dk 3 d k4 ■ 

This corresponds to taking only the on-site interaction in 
the d-orbital into account. Since U' = 0, the self-energy 
insertion contributions to all four running couplings take 
on the same value and the direct particle-hole diagrams 
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FIG. 6. (Color online) Flow in the low-energy model for AE = 19.5, w d = 11.5, w s = 4.5, U = 0.03 • -liv 2 w d , U' = (a) and 
U' = 0.317 (b) at approximation level 1. The two-patch couplings are depicted as solid (gi), dashed ($2), dashed-dotted (33) 
and dotted lines (54). The band-flip self-energy term in the feedback has not been taken into account. 
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FIG. 7. (Color online) Flow without band-flip self-energy 
term for AE = 6.5, w d = 11.38, w s = 23.51, U = 0.03-47r 2 w d , 
U' — 0.3(7 at approximation level 1. Line styles as in Fig. |S] 

(IB4[) vanish. In the diagrams in Eq. (IB2|) the hybridiza- 
tion Cfc appears only inside the integrand whereas the 
external legs of momentum k have a factor dk, which 
in contrast to Cfc is invariant under a spatial rotation 
by 7r/2. This implies that the particle-particle contri- 
butions to AV-(A,A, B) and AV- (A, A, A) are of equal 
size, whereas the crossed particle-hole diagrams give dif- 
ferent contributions. The latter can be seen as follows: 
After calculating the Matsubara sum for A = 0, the inte- 
gral in Eq. (|B3I) reads as 

/ dqcf+X [S+fa+lJ-iMq)] -1 © (-£-(<!)) > (8) 
Jbz 

with the Heaviside function Q(x). For AVL(A,A,A), 
which renormalizes 34, we have 1 = and for 
AV—(A,A,B), which renormalizes 53, the transfer mo- 
mentum is 1 = 7r. Since both bands have a curvature with 
the same sign in our model and since E + is always posi- 




Qt 



FIG. 8. (Color online) Phase diagram of the t-t' Hubbard 
model in two-patch approximation, approximation level 1, 
with qt = —t'/t. Points that correspond to the same value 
of U e s are connected by a line. The region around the sep- 
aratrix between the dSC and FM instabilities is inaccessible 
in our approach, since the numerical integration of the loops 
gets too cumbersome. 



tive, the denominator in the integrand of Eq. (|5J) should 
take on smaller values for 1 = than for 1 = # on a 
large phase space region centered around q = 0. One 
might therefore expect a suppression of (73 — 34. This ar- 
gument, however, ignores the momentum dependence of 
the orbital weight completely. Whilst being zero along 
the diagonals of the BZ, the hybridization matrix ele- 
ments c q have their maximal value close to the saddle 
points of the lower band. This weakens the effect of the 
band curvature on 33 — 34 , in particular if the hybridiza- 
tion shows plateau-like structures centered around the 
van Hove points as for the parameters underlying the 
flow in Fig. [9] At these points, however, the four-fold ro- 
tation symmetry of the dispersion gives rise to identical 
denominators of the loop integrand in Eq. © for 1 = 
and 1 = iv. Hence, the effect of band curvatures that 
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one might expect ignoring orbital makeup effects should 
play a minor role in such a case. In contrast, the phase 
space weight imposed by the orbital makeup may lead to 
an enhancement of 53 — (74 , even for weak hybridization 
as in the case of Fig. [9j For larger values of the hy- 
bridization, the situation is yet more involved, as terms 
with opposite signs compete. So far, we have discussed 
the contributions of AVI to 173 — 174 at approximation 
level 3. We have found that they are quite sensitive to 
the orbital makeup. The question now is how they affect 
the critical scale: On average, AV_ suppresses the cou- 
plings gi while the difference g% — #4 may be enhanced 
or lowered depending on the orbital makeup. In general, 
smaller initial values of the couplings give rise to lower 
critical scales while a larger value of the d-wave coupling 
.93 — .94 promotes a sign change of 174 at higher scales. 
So for enhanced d-wave coupling, we have to deal with 
two counteracting tendencies and it is not clear a priori 
which one prevails, whereas a suppression of A cr ;t is to be 
expected if 33 — 174 is lowered. In Fig. [TUIthev lead to a 
surprisingly large enhancement of A cr it ■ Indeed, (73 — 54 
gets larger when the six-point feedback is taken into ac- 
count, but if we neglect contributions from d,2 and d% in 
the flow equations Eqs. (|4][7]) the critical scale is virtually 
the same as in the conventional truncation. Moreover, its 
value is increased by orders of magnitude without the di 
and g?3 terms. This points out the importance of those 
terms and their interplay with the six-point feedback for 
the dispersion underlying Fig. |H1 

In Fig. [10] 53 — #4 is suppressed by AV_. Without 
the d,2 and c?3 terms, this would only change the critical 
scale by values below the level of accuracy, while A cr ;t 
is significantly enhanced if those terms are taken into 
account. As for the flow to a ferromagnetic instability 
in Fig. we find that the six-point term is of minor 
importance for the corresponding parameters. 

Finally, we investigate the impact of the band sepa- 
ration AE on the critical scale as depicted in Fig. ITT1 
Lower values of AE correspond to a larger ratio —t'/t and 
therefore to a lower critical scale in the dSC regime. We 
observe that for band separations smaller than 19.5 units, 
the six-point feedback substantially enhances the critical 
scale. In particular, when the critical scale gets small due 
to the competition with the FM channel, the six-point 
term can change the result by an order of magnitude. 
This indicates that these corrections may play a role in 
situations with competing ordering tendencies. Deep in 
the d-wave regime, or also on the FM side, the impact of 
the six-point term is only of quantitative nature. Since 
the behavior at small dSC critical scales mainly stems 
from the interplay of the six-point feedback and the di 
loop, it should result from the increasing contribution of 
c?2 that we have to encounter when AE is lowered. The 
inclusion of the band-flip self-energy term only leads to 
slight changes of the critical scale and does not affect our 
result on a qualitative level. 



VI. CONCLUSIONS 

In summary, we have proposed an fRG scheme for 
multiband systems that takes virtual excitations outside 
a low-energy window into account. This approach in a 
way extends the usual RG approach of accessing the low- 
energy physics of a given system, where the effective ac- 
tion is truncated after the four-point interaction. The 
conventional treatment completely neglects virtual exci- 
tations involving both high- and low-energy modes once 
one has switched to the low-energy theory. Our approach, 
in contrast, considers them up to second order in pertur- 
bation theory within a frame with a larger energy scale 
around such a low-energy window. On the formal level, 
this is achieved by a one-loop fRG treatment of an effec- 
tive action for the low-energy modes, which is truncated 
after (i.e. keeps) the three-particle interactions generated 
by the high-energy modes (six-point term). Of course the 
low-energy theory could also be solved using a different 
method than fRG, and also here the six-point term of the 
effective interactions may play a role. 

In the RG flow in the low-energy window, the six-point 
vertex gives rise to mixed one-loop diagrams with one 
low-energy and one high-energy leg. Actually, a part of 
these diagrams is summed up in the cRPA framework 
currently used in ab-initio calculations^. The correla- 
tion functions of high-energy modes are no longer kept 
track of in our scheme. This drastically reduces the re- 
quired numerical resources compared to a 'full' fRG ap- 
proach with an extended energy window. 

We have numerically investigated the impact of the 
mixed diagrams that are absent in the conventional ap- 
proach for a simple two-band model within the two-patch 
approximation. Leaving a large part of the Fermi sur- 
face unaccounted, this approximation is mainly chosen 
for reasons of technical simplicity. On the other hand, 
the two-patch treatment may appear questionable if vir- 
tual excitations in the upper band away from the Fermi 
level are considered, since renormalization effects at in- 
termediate scales are discarded. However, the impact of 
the rest of the Fermi surface has been investigated for the 
one-band Hubbard model within iV-patch schemes 10 , and 
recently within a channel decomposition 2 -. The leading 
tendencies in the parameter region of interest are well 
captured in the two-patch treatment. So our results may 
indicate well how the virtual high-energy excitation affect 
the RG flow, also if the full Fermi surface was considered. 
A thorough fRG study of multiband models that takes 
into account the off-Fermi-surface bands perturbatively 
would require a multi-patch scheme that goes beyond 
Fermi-surface projection of the interactions. The authors 
are currently implementing such a scheme. 

In this work, we have used a multiplicative frequency 
cutoff 20 , which allows for a second dSC and a FM fixed 
point that cannot be found in two-patch studies using 
a traditional momentum-shell cutoff. This is due to the 
particle-hole diagrams with vanishing momentum trans- 
fer that contribute to the flow at nonzero scales if a mul- 
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FIG. 9. (Color online) Flows in the low-energy model for AE = 19.5, w d = 11.5, w s = 4.5, U = 0.03 • A-n 2 w d , V = at the 
different approximation levels. The interaction has been projected to the 52/173-54/53 plane, while the ^-direction corresponds 
to the scale. The impulse-type vertical lines have been added for clearness and have no physical meaning. The initial flow for 
high scales is shown on the left hand side. Here one can see that approximation level 2 with six-point feedback has the same 
initial condition for the 5;s as level 1, but approaches quickly level 3 in the early flow. On the right hand side, the continuation 
of the flow at lower scales is depicted. Level 2 and 3 are basically equivalent and have a higher critical scale than level 1. 
The ratio 52/53 flows to zero, corresponding to the second d-wave pairing fixed point discussed in Subsec. [V Al The band- flip 
self-energy term in the feedback has not been taken into account. 

level 1 level 1 




FIG. 10. (Color online) Same as in Fig. H but for U' = 0.317. 
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AE 

FIG. 11. (Color online) Critical scale A cr it for Wd — 11.5, 
w a = 4.5, U = O.OS-lvr 2 ^, U' = 0.3(7 as a function of AE, in 
the dSC regime. Circles correspond to approximation level 1, 
diamonds to level 2 and crosses to level 3. Green and black 
markers represent values obtained with and without band-flip 
self-energies, respectively. 



tiplicative frequency cutoff is used. The regularization 
scheme of momentum-shell type put forward in Ref. [26| 
also accounts for these diagrams and gives results that 
are qualitatively similar to ours at approximation level 1. 
However, it cannot be extended to a region outside the 
patches and is therefore not viable at approximation 
level 2. 

In the two-patch treatment the approximation levels 
we call 2 and 3 merely coincide, that is the six-point 
term can be described by adding the mixed diagrams to 
the bare values of the couplings in the initial conditions 
of the flow. This represents a technically appealing ap- 
proximation that can be implemented more easily com- 
pared to keeping a flowing six-point term. It should hold 
in cases where the mixed diagrams are renormalized at 
scales at which the low-energy diagrams only flow weakly. 
Close to a separatrix in the flow between dSC and FM 
ordering, even slight changes in the parameters of the 
bare action, i.e. the initial conditions of the flow, can re- 
sult in a drastic change of the critical scale. In the two- 
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patch model, this can be understood by the small number 
of fixed points with distinct physical properties, which 
can be mutually exclusive. In an TV-patch approach this 
should in principle no longer hold, since there could be 
a larger number of fixed points that decay into (partly 
overlapping) classes sharing similar physical properties. 
Yet, the flows for the full Fermi surface still clearly show 
a close competition and a possible quantum critical point 
between dSG and FM instabilities. Indeed, the sensitiv- 
ity with respect to the initial conditions is strongly re- 
flected in our study of different approximation levels to 
the effective theory. We find that the six-point feedback 
can enhance the critical scale by orders of magnitude in 
the dSC regime, provided that the particle-hole bubble 
with zero transfer momentum diverges strongly enough 
in the infrared, i.e. there is a strong competition between 
dSC and FM tendencies. Away from this critical region 
and on the FM side, the impact of the six-point term is 
reduced. 

Currently it is unclear if these differences will persist 
or get weakened if a multi-patch scheme is used, and if 
the six-point feedback may still play an important role at 
least on a quantitative level in such a flow. Very recent 
results on two- and three-band models aiming at different 
trends in cuprate systems using related approximations 
gave a slight enhancement up to 15% due to the bands 
away from the Fermi leve£. Regarding the iron pnictide 
superconductors, there a several very interesting ques- 
tion that arise from this study. There it might be useful 
to carefully study other observables than just the critical 
scale. For example, most RPA and fRG studies indi- 
cate that the orbital character of the bands induces a 
pronounced anisotropy of the sign-changing s-wave gap 
function around the various Fermi pockets. Of course, 
for the predictive power of these calculations it is neces- 
sary to check how these anisotropics are affected by the 
higher-order terms in the effective interactions. Then, 
depending on the parameters of the multiband model, 
there is the possibility of a competition of other supercon- 
ducting states with the anisotropic s-wave state. Again, 
in such situations additional effects may play a decisive 
role, and may either increase or decrease the degree of 
competition. 

In summary, we have provided tractable approxima- 
tions for improved effective interactions within the con- 
duction band of a multiband problem and performed a 



case study in which the correction terms play a visible 
role. The next steps should include more extended appli- 
cations in multiband models for correlated electron sys- 
tems, in order to understand the importance of these 
corrections in a comprehensive way. 
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Appendix A: Parameterization of the interaction 

We will now parameterize the interaction in a way sim- 
ilar to±l The quartic part of the action is therefore cast 
in a suitable form: 

Ct,{ki ,CTj} 

U' 

+ — ^ <5qi,-q 2 [E a E a — Da-D a ) 8k 

{ai,ki ,cr;} 

where D a = S aua4 S a2 ^ 3 , E a = 8 ai>(T3 5 a2!ai , and where 
5f. = S (ki + k- 2 — fc 3 — k^) ensures energy and momen- 
tum conservation. 

Upon the orbital-to-band transformation the 5- 
function in the momenta gets multiplied by a product 
of form factors dk and/or Cfc. We now rewrite S^ 4 -* in 
terms of the fields \ 

S(i) ix,x}=\j2 F(Xi,X 2 ,X 3 ,X 4 ) 

{Xi} 

x x(x 1 )x(x 2 )x(x 3 )x(x i ) 

and decompose the interaction according to the band in- 
dices of the external legs. This yields 



F(Xi,X2.,X 3 ,X i ) =Ax 1 ,X 2 Ax 3 ,X 4 . [^+(6:^2,^3,^4) ^a, ++++ + ^-(^1,^2,^3,^4) 

+ F + ^-fa,£ 2 ,&,U)8 a ,-- + + + F-^fa,£2,&,&)6 a , + +-- +4 J F±^±(a,6,6,£4)<W— + 

+ 2F + ^ ± (a,6 1 6,C4)^,-++ + + 2 J F_^ ± (a,6,6 ) ^)5«,+— 

+2F±_> + (&, *«*,+++- + 2F±^_(£ 1 , £3, U)5 a< +] , 



I 

where £j denotes {ki,(Ji) and Xi = (cki,£i) and where the two-point antisymmetrization operator A has been de- 
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fined as A a ,bf{a, b) = [f(a, b) — f(b, a)] /2. In the orbital 
picture, the interaction is invariant under a band-index 
flip. Since the trace of the matrix of the orbital-to-band 
transformation Eq. ((5J) vanishes, this property also holds 
in the band language, giving rise to the following identi- 
ties 



*±->±(6, £ 2 ,£ 3 ,£ 4 ) = F ± _ ± (6, £l,£4,£ 3 ) 
>±(£l,6,6»,&) = ^±(6, 

F±^-(6 J 6,6^4) = i ? ±^+(?i,6,6,C4). 

Hermiticity of the Hamiltonian in the band language re- 
quires that the band conserving terms F_, i<±_>± 
must obey 

6, = (Al) 



For the non-band-conserving terms it leads to the follow- 
ing relations: 



reflected by the following antisymmetry constraint 



= F:^ + (^,e 3 ,6,a) 



The anticommuting nature of the Grassmann fields is 
I 



-Fx(£l,£2,£4,£ 3 ) , 



for X = +, — , 2. According tois, these U(l)-vertices can 
be parameterized as 

=6k [E rT V x {k 2 MM 
-D <T Vx(k 1 ,k 2 ,k 3 )} , 

where Vx(ki, k 2 , k 3 ) = Vx{k 2 , ki,k 3 — fci — /c 2 ). For ver- 
tices with three legs on one band and one on the other, 
however, one of the symmetry constraints is violated and 
we only have F s (£ x , £2 , £3 , £4) = ~F 3 (£1 , 6 , £4 , £3) , which 
allows as well for the parameterization 

i ? 3(£i J £ 2 ,£3,£4) =S k [E^ktMM) 
-ArVaC&i.fc,**)] • 

In contrast to the vertices Vx considered before, there 
is no symmetry constraint for V 3 . Finally, there are no 
antisymmetry relations for which gives rise to the 

parameterization 



f±^±(£i,£2,£3,£4 



with the symmetry constraint reflecting the behav- 
ior under a band-index flip. V±^±)(ki,k 2 ,k 3 ) — 

V±^l± (k 2 ,ki,k 3 — ki — k 2 ). We now give explicit ex- 
pressions for the coupling functions 



y_ (fci , fc 2 , fc 3 ) = u (j[ d ki + n Ck ^j + u' (d kt , 



v c k 2 c k 3 d ki + c kl d k2 d k3 c k4 ) 
V 2 (ki,k 2 ,k 3 ) = U (d kl d k2 c k3 c ki + c kl c k2 d k3 d k J - U' (d kl c k2 d k3 c ki + c kl d k2 c k3 d ki ) 
Vl°} ± (ki,k 2 ,k 3 ) = U (d kl c k2 c k3 d ki + c kl d k2 d k3 c k4 ) + U' ^JJ dfc, + JJc fe< J 

V^l ± {k u k 2 ,k 3 ) = U {d kl c k2 c k3 d kl + c kl d k2 d k . 3 c kl ) - U' (c kl c k2 d ka d ki + d kl d k2 c k . 3 c k J 
V 3 (k 1 ,k 2 ,k 3 ) = U (d kl c k2 c k3 c ki - c kl d k2 d k3 d ki ) + U' (d kl d k2 d ks c ki - c kl c k2 c ks d ki ) . 



(A2) 



(A3) 



So the bare interaction can be expressed in terms of 
5 independent functions of three 1 + 2-momenta by ex- 
ploiting its symmetries. 



Appendix B: Flow equations 

In this appendix, the flow equations used in this paper 
shall be given explicitly. We neglect self-energy effects 
and put the six-point vertex to be constant during the 
flow. Therefore, only the four-point flow equation mat- 
ters. If the six-point function is ignored, the right hand 
side of this equation is given by Eq. (88) in Ref. [13 with 
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V- as an initial condition for V. In this form, the flow 
equation can be used to discuss the impact of orbital 
makeup. Proceeding further, we take into account the 
feedback term 

A 4 (a,6,6,a) = -\f dnS{m,m) 

xF (6) (»h,Ci,6, 
A 4 (6 6k [E* 5V- (hi MM) 



S denoting the single scale propagator S = G — GSG 
with self-energy S. In the following, we assume the cou- 
pling functions V_ and V 3 to be real as is the case for the 
model analyzed in this work. Up to second order in the 
bare interaction, the six-point coupling function is given 
by 

F^iti ...&)= i; ..; : -4; ; .... J d v G+( m , m ) 

xf8(iii,6,a&)% 1 6 ) 6 1 &) , (Bi) 

where G + denotes the upper-band propagator and where 
the three-point antisymmetrization operator 



A a ,b.c.f{a,b,c) = — ^P w /(7r(o),7r(6),7r(c)) . 



is given by the difference of the sums of cyclic (P n = 1) 
and anti-cyclic (P w = —1) permutations n. These anti- 
symmetrization operators give rise to self-energy inser- 
tion and one- loop diagrams A 4 = A SE + Ai oop . The 



former contain the band-flip self-energy 

£±(£1,62) = <5(£i-&) JdnidrfrF z (£ u vi,rfr,&)S(vum) 
and read as 

A S e(£i,6,£3,£4) = Jdi] 3 dr] A G + {-q 3l -q A ) 

X i ,£a F 3 (% , & , 6 , £4 ) S± (»?4 , £l ) 

+ ^?3,«4^3(??4,6,6,£l)S±(7?3,£4)] ■ 

This expression can again be parameterized as 

A SE (£i, =4 [S CT ^ E (fc2,fci,fc 3 ) 
-^ CT ^SE(fcl,fc 2 ,fc 3 )] , 

where Vse obeys the same symmetry constraints as V_. 
This gives rise to 

VsE(fcl,fc2, fc 3 ) = V r 3(fcl,fc2,fe3)S±(fcl) 

+ y 3 (fc 2 ,fc 1 ,/c 1 + /c 2 -fc 3 )s ± (/c 2 ) 

+ V 3 (k 1 + k 2 - k 3 , k 3 , fc 2 )S±(fci + fc 2 - fcs) 

+ y 3 (fc 3 , fci + fc 2 - h , fci )s± (fc 3 ) 

with 

£±(k) = G+(k) J dqS(q) [V 3 (k,q,k)-2V 3 (k,q,q)] 

The one-loop part comprises particle-particle, crossed 
and direct particle-hole diagrams 

Aioop (£1 , 6 , £3 , U) = h [Ea- Vioop (k 2 ,k 1 ,k 3 ) 

-Da- Vl 00 p(fci,fc 2 ,/C 3 )] 
Vloop(fcl, &2, A3) =^pp(fcl,fc 2 ,^3) +7?-phxr(fcl,fc2,A: 3 ) 

+ ^ p h,d(fci,fc2,A: 3 ). (B2) 

Since V3 obeys no symmetry constraint, the particle- 
particle contribution 



n pp (h,k 2 , k 3 ) = -J dq S(q) G+(l - q) [V 3 (l - q, q, h)V 3 (l -q,q,l- k 3 ) 

+ V 3 (l- q, q, k 2 )V 3 (l - q, q, fc 3 )] /=fel+/£2 

consists of two terms that otherwise would coincide. The crossed particle-hole terms 

H ph , CI (k 1 ,k 2 ,k 3 )=- JdqS(q) G+(l + q) V 3 (l + q,h,q)V 3 (l + q,k 2 - Z,«)| I=fc ,_ fcl 

- JdqS(q)G + (l + q)V 3 (l + q,k 2 ,q)V 3 (l + q,k 3 ,q)\ l=ki _ k3 (B3) 
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behave likewise. The direct particle-hole diagrams read as 



Note that for vanishing inter-orbital interaction U' = 
the direct particle-hole diagrams vanish due to 
V 3 {k 1} k 2 , fc 3 ) = V 3 (ki,k 2 , k 1 + k 2 ~- k 3 ). 



Appendix C: Two-loop corrections 

Here, we give an argument why not going beyond 
level 2, i.e. why dropping the renormalization of the six- 
point vertex itself, may suffice. If one had not done so, a 
correction term would have to be added to the constant 
three-particle vertex as in Fig. IT2"a). By integrating and 
then iterating the flow equation for the six-point vertex, 
this renormalization induced term can be expressed as 
a sum of diagrams containing the four-point and con- 
stant six-point vertices only up to arbitrary order. Note 
that the constant part of the six-point vertex is of sec- 
ond order in the four-point couplings. In leading, i.e. 
third order, we obtain two one-loop terms depicted in 
Fig. 112b). One of these diagrams includes only four-point 
vertices and would as well be present in the absence of 
an initial six-point term while the other one contains this 
initial three-particle interaction. When the right-hand 
side of Fig. [Fib ) is now fed back into the flow equation 
for the four-point vertex, diagrams with overlapping and 
non-overlapping loops arise. If we had started our fRG 
analysis directly from the full action (S^ + S^ 4 * 1 ) [x, x] 
instead of the effective low-energy action S e g and kept 
the band indices as variables attached to the legs of the 
vertices, the overlapping ones would be neglected in the 
Katanin truncation. So they shall be dropped in the 
present treatment as well. 

We then end up with the correction terms depicted 
in Fig. 1131 The first and the third one can be merged 
with the one-loop terms in Fig. @] leading to a Katanin 
substitution S — > G both in the low-energy term as in the 
one-loop feedback. Since we neglect self-energy effects in 
our numerics, this substitution will not change our results 
for the six-point feedback. 

The second diagram in Fig. 1131 however, requires more 
care. We now proceed with giving upper estimates for the 
remaining correction term and the one-loop feedback and 
low-energy terms in Fig. |4j If the frequency dependence 
of the coupling functions is dropped, all Matsubara sums 
can be evaluated analytically giving rise to the following 
rules for an estimate. 



(B4) 

I 




b) 

FIG. 12. a) Full (scale dependent) vertex (hatched circle) 
written as a sum of its initial value (empty circle) and a 
renormalization induced correction term (hatched square), b) 
Leading order result for the correction term. 




FIG. 13. Leading corrections to flow equation in Fig. [4] Dia- 
grams with overlapping loops have been neglected. 

• The four-point coupling functions V 3 and V- at 
scale A are replaced by their maximal value g 3 and 
<?_, respectively. 

• Mixed loops including a scale derivative are re- 
placed by a factor 47r 2 (e+A) _1 , where e+ denotes 
the minimal energy of the high-energy bands. The 
band-flip self-energy diagrams behave likewise. 

• Lower-band loops with and without a scale deriva- 
tive are replaced by 47r 2 A~ 2 and 47r 2 A~ 1 , respec- 
tively. 

At all scales the correction term Fig. [T3] should be small 



ft P h,d(fci, k 2 , k 3 ) = - J dqS{q) G+(l + q) [-2V 3 (l + q, h,ki + l)V 3 {l + q, ha, k 2 ) 
+ V 3 (l + q, k u ki + l)V 3 {l + q, k 3 , q) + V 3 {1 + q, h,q)V 3 (l + q, fa,k 2 )] l=k2 _ 
-Jdq S(q) G+(l + q) [-2V 3 (l + q, k 2 , k 3 )V 3 (l + q, fc x - I, h) 

+ V 3 (l + q, k 2 ,ki - l)V 3 (l + q, k 3 , q) + V 3 {1 + q, k 2 , q)V 3 {l + q,h~ I, ki)] h 
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compared to the the low-energy term, which implies 

4tt 2 ^ « <?_£+ . (CI) 

Note that the orbital makeup reduces the value of the 
bare coupling functions and may therefore finally allow 
for the omission of the correction term. At scales at 
which the one-loop feedback term flows, the one-loop 
feedback should prevail against two-loop corrections, i.e. 
47r 2 g_ <C A. Together with the condition Eq. (jCip , this 
requires the cutoff A to be much larger than <?§/e + . At 
lower scales A <C <? 2 e + /g|, however, the mixed one-loop 
diagrams eventually saturate and the feedback term be- 
comes negligible compared to the low-energy loop term 
that may finally drive the flow to a strong coupling fixed 
point. The crossover region between these two regimes 
should be small, as long as the inequality (jClj) holds. So 
a one-loop fRG approach should suffice to qualitatively 
discuss the impact of the six-point term on the critical 
scale for d-wave superconductivity, for example. 

Finally, we feel that a comment on the relation of the 



flows of r and T_ is in order. If self-energy effects are 
neglected completely, the four-point vertex of T_ is equal 
to the four-point vertex of T with all external legs on 
the low-energy bands, since they both must lead to the 
same four-point correlation functions for the low-energy 
modes. If one now considers the flow of T for an energy 
shell cutoff in the usual truncation and in addition forces 
all four-point vertices with at least one leg on the high- 
energy bands not to flow, the one-loop feedback in the 
flow of T„ is recovered. In the RG flow of F, the most 
relevant correction term to this approximation consists 
of the diagram with the internal lines on the low-energy 
bands and three external legs on the low and one on the 
high-energy bands. If this correction is fed back into 
the flow of the vertex with all external legs on the lower 
bands, the feedback correction term is equivalent to the 
second diagram in Fig. [13] in leading order. As soon as 
self-energy effects are taken into account, however, this 
correspondence breaks down. 
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